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We derive the modes inside a cylindrical waveguide of finite surface 
impedances, assuming the waveguide transverse dimensions are 
large compared to the wavelength \. This paper restricts its consid- 
eration to the modes with ft ~ k, where ft is the propagation constant 
and k = 2-rr/X. For these modes we show that asymptotically, for large 
values of k, the field \f/ becomes infinitesimal (of the same order of 
1/k) at the boundary. Taking this into account, we obtain simple 
expressions for the asymptotic properties of $ for large k. The theory 
applies to a variety of waveguides: corrugated waveguides, optical 
fibers, waveguides with smooth walls of lossy metal, and so on. An 
important property of the modes considered here is that their atten- 
uation constant is very low, i.e., asymptotic to 1/k 2 for large k. Thus, 
these modes are useful for long-distance communication. Another 
consequence of $ — > at the boundary is that for large k the 
distribution of$ inside the boundary is essentially independent of the 
boundary parameters, i.e., independent of the surface impedances in 
the longitudinal and transverse directions. This consequence implies 
that the same radiation characteristics of the corrugated feed can be 
obtained using other structures and, therefore, construction can be 
simplified in many cases, with little sacrifice in performance. We also 
derive general expressions for ty and the propagation constant ft. 

I. INTRODUCTION 

It is known 1 " 8 that in certain waveguides the field becomes, under 
certain conditions, very small at the boundary. Consider, for instance, 
a corrugated waveguide of radius a and let X be the free-space wave- 
length. This waveguide is characterized at the boundary by finite 
surface impedance Z z in the longitudinal direction. The frequency 
dependence of Z 2 , which is determined by the depth of the corruga- 
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tions, causes the transverse field distribution \p(x, y) of a mode to vary 
with the frequency k = 2tt/\. However, this frequency dependence 
virtually disappears (for all modes except surface waves) if the wave- 
guide dimensions are large enough. In fact, one finds that ${x, v) 
approaches for ka -* oo a frequency independent distribution that 
vanishes at the boundary. 5 This behavior is responsible for the low 
attenuation constant, for the excellent radiation characteristics, and 
the wide bandwidth of corrugated waveguides. 5 We show here that the 
same behavior also occurs, under quite general conditions, in a variety 
of uncorrugated waveguides. 1 " 21 Figures la and lb show two examples, 
a dielectric waveguide 7 ' 813 " 16 of general cross section and a hollow 
waveguide with metal walls coated by a dielectric layer. 417 Other 
examples can be obtained by modifying the boundary conditions in a 
variety of different ways. For instance, several dielectric layers may be 
used in Fig. lb, or a metal grid of transverse wires may be placed at 
the boundary, as pointed out in Section II. Other examples are the 
waveguides of dielectric or lossy metal considered in Ref. 2. We now 
outline the main assumptions. 

Consider a cylindrical waveguide with an interior region of homo- 
geneous and isotopic material, as in Fig. lc. Let Z and k be the wave 
impedance and propagation constant for a plane wave in the interior 
region, and let C denote the boundary. Then 



-4 



k = coveu. (1) 



Consider a mode with propagation constant fi and let E, denote the 
transverse component of the electric field, 

Ei = «K*, y)e~ jfi '. (2) 

Let 2a be a characteristic dimension of the waveguide, for instance the 
width in the x direction as in Fig. lc. We are concerned about the 
asymptotic behavior of \p(x, y) for large values of ka. Consideration 
will be restricted to the modes for which the propagation constant /? 
approaches k, as ka — » oo. Thus, we assume 

/? -► k for ka -> oo. (3) 

This excludes surface waves, as pointed out in the following section. 
Then, a property of the modes considered here is that E and H become 
transverse, in the limit as ka — * oo, 

lim E z - H z = 0. (4) 

to >co 

Another property is that asymptotically, per large ka, a set of linear 
relations exist at the boundary among the tangential components of E 
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TANGENT PLANE 



Fig. 1 — Three examples of cylindrical waveguides: (a) a dielectric rod, (b) a waveguide 
with metal walls coated by a thin dielectric layer, and (c) a waveguide with boundary 
conditions shown in Fig. 2. 



and H. It is convenient to write these relations in the form 



ZH T 



-J[H] 



-ZH Z 



on C, 



(5) 



where [H] is a 2 x 2 matrix and H T , E T denote the components of H, 
E in the direction of the unit vector r in Fig. lc. These relations 
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together with eq. (4) give at the boundary 

lim E r = H T = on C, (6) 

te-MO 

provided the matrix [H] does not diverge for ka — > oo, 

[H]*oo for ka -»». (7) 

Throughout this article, we assume conditions (3), (5), and (7). 
Condition (5) is discussed in the following section, where it is pointed 
out that for most waveguides considered here, [H ] is a diagonal matrix. 
In this case it is convenient to define at the boundary surface imped- 
ances Z z and Z T by writing 

E r = Z T Hz, E z = —ZzH T , (8) 

where Z z is the longitudinal impedance, and Z T is the transverse 
impedance. Then condition (7) requires 

Z T , l/Zz * oo for ka -> oo. (9) 

It is important to realize that this requirement is violated in a number 
of cases. It is violated in a hollow waveguide with metal walls of perfect 
conductivity, since then Z 2 = 0. Furthermore, in a corrugated wave- 
guide with corrugations of depth d, the longitudinal impedance Z z is 
determined by kd, and there are certain frequencies for which Z z — 0. 
A similar situation arises in Fig. lb where both Z T and Z z vary with kd. 
Throughout this article it will be assumed that the quantities 

1 U J_ Z_ 

ka' ka Z' kaZz 

are small. Therefore, the results will not apply in the vicinity of the 
above frequencies. 

A direct consequence of condition (6) is that the boundary values of 
\f/ vanish, in the limit as ka —* oo, 

lim 4>(x, y) = on C. (10) 

ka-icc 

Another consequence is that \\i approaches a distribution i//*, indepen- 
dent of ka, for ka —* oo. For finite ka, 

+ = +~ + 8+, (11) 

where 5»// (but not $» ) varies with ka and 

lim $V - 0. (12) 

Ad-MD 

Notice condition (10) implies that 

4*. = on C. (13) 
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These results are of practical interest for several reasons. In the design 
of a feed, 6 it is desirable that the boundary values of $ be small [as 
implied by eq. (10)] since these values determine radiation in the side- 
lobes due to edge diffraction at the aperture. Furthermore, for broad- 
band applications, it is desirable that the frequency dependent part 
&// of the aperture illumination be small, as implied by eq. (12). Finally, 
in a corrugated waveguide, or a waveguide with metal walls coated by 
dielectric layers, power is lost only at the walls and, therefore, it is 
determined by the boundary values of \p. Since these boundary values 
vanish for ka — > oo, the attenuation constant for the above waveguides 
for large ka is very small 1,2,418 " 21 ; it is asymptotic to ka~ 2 . We shall see 
that in general *//« is independent of [H] and, therefore, a variety of 
different waveguides, with different [H] but the same boundary shape, 
give rise to the same \f/„ . This explains the similarity, noted in Ref. 9, 
between the modes of a corrugated waveguide and those of an optical 
fiber, a dielectric lined waveguide, 4 and a hollow dielectric waveguide. 2 
This similarity implies that essentially the same radiation character- 
istics of corrugated waveguides can also be obtained with a variety of 
other waveguides. 

In the particular case of the optical fiber, some of our results are 
implied by the asymptotic expressions derived in Ref. 7. Exact solu- 
tions for the modes of the corrugated waveguide, 12 the optical fiber, 7,8 
and the hollow waveguide of dielectric 2 are known for circular geom- 
etry. For a rectangular cross section, only approximate solutions 3,22 are 
known, except in special cases. 23 Exact solutions for the slab waveguide 
are given in Refs. 8 and 24. In all these cases one finds that condition 
(3) implies condition (10). Measurements of the radiation characteris- 
tics of a dielectric horn are described in Refs. 25 and 26. Some of 
the properties derived here apply also to propagation in stratified 
media. 27 " 29 The use of surface impedances to characterize a boundary 
is discussed in Ref. 30. 

II. BOUNDARY CONDITIONS FOR ka -> oo 

We now derive and discuss eq. (5). Figure lc shows a waveguide 
directed along the z axis and of general cross section in which v is the 
outwardly directed normal and t is a unit vector tangential to the 
boundary, 

t = i, X v. (14) 

The medium inside the boundary is assumed to be homogeneous and 
isotropic. Let C denote the closed contour of the boundary in the plane 

2 = 0. 

We are concerned with the properties of \}/{x, y) in a waveguide of 
large transverse dimensions. Thus, consider a mode propagating in the 
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waveguide of Fig. lb and suppose the width 2a is increased keeping 
the dielectric thickness d fixed. The resulting behavior of $(x, y) as ka 
— > oo can be derived exactly in two cases, when (see Appendix C) 

^ = (15) 

dy 

and when the boundary is a circle. In both cases one finds that for 
most of the modes /8 — ► k, as ka -* oo. For these modes, the normalized 
field amplitude \p(x, y)/*p(0, 0) becomes infinitesimal at the boundary 
for ka — » oo. For the remaining modes, those for which /3 does not 
approach k, just the opposite behavior takes place: The field becomes 
confined to the immediate vicinity of the walls, degenerating into a 
surface wave with propagation constant determined by the surface 
impedances of the walls. Here, consideration will be restricted to the 
modes satisfying condition (3). An important property of these modes 
is that asymptotically, for large ka, the surface impedances Z r and Z z 
become independent of ka. In fact, if one writes 

*-4 y ~4 < i6 > 

then in Fig. lb 



y^ tan(vV^lAd), (17) 

Y->-^£= . } (18) 

v^^T tan(V/i 2 - 1 kd) 

as shown in Appendix C. Thus, Z T and Z z depend only on the refractive 
index n and the thickness kd of the dielectric layer. 

For a circular boundary, the above relations can be derived rigor- 
ously by expressing the field in terms of Bessel functions, and then 
making use of well-known expressions giving the asymptotic behavior 
of these functions for large arguments as in Ref. 5. They can also be 
derived by the following argument, which applies in general to a 
boundary of arbitrary shape (Fig. lc). Consider the field in the vicinity 
of a boundary point P in Fig. lc. Since ka is large, the curved boundary 
can be approximated locally by the tangent plane. Furthermore, since 
/? = k, the field is produced locally by a plane wave at grazing 
incidence. If one determines the law of reflection for such a plane 
wave, and takes into account that the plane of incidence is parallel to 
the z axis, one obtains 24,27 eqs. (16) -(18). 

The above argument can be used to derive the asymptotic behavior 
of Z z and Z T for a variety of other waveguides, illustrated in Fig. 2.* In 
case (d) the boundary is a smooth surface of lossy metal with surface 



* Notice ko denotes k in free space. . 
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resistance R s . In (a) the metal surface is corrugated. In (b) the medium 
outside the boundary is dielectric with refractive index n-i < rii. This 
corresponds to the optical fiber of Fig. la. In (c), rc 2 > »i as in Ref. 2, 
and, therefore, both Z z and Z T are real, which implies power is lost 
through the boundary. Other boundaries of practical interest are 
obtained from Fig. 2 by placing a grid of thin wires tangent to t on the 
boundary. This will cause Z T — in all cases. For all these waveguides 
one finds, for large ka, that [H] in eq. (5) is a diagonal matrix, 



[H] = - 



Y 
X 



(19) 



TYPE OF BOUNDARY 



-■*— V 



T 




rii 



jz(l- JLjTANK d 
( h « X ) 



Zt-0 



(b) 



(ni > n 2 



v< 



?-i 



n 2 



(z- -\JjI77) 






(c) 




(r = Vn|-n?/n 2 



n,< n 2 » 



1 "1 




(i+i)R s 



(1+JJR, 



I METAL 

-I V,i,'(/,',, , l<U . f<iW<<> ' 



i'lrr ZT + R.O+T' 




(T = TANK Vn 2 -n? d, 
-^r- « T. 1/T) 



■ ' "1 



+ R.H + T 2 



(e) 



(n 2 > n,) 




'^- Z 



T' + |T 
1 + jTT' 



_L JlL z T + jT 



r n 2 1 + jTT' 



V n 3 - n 1 ^2_ 



v^ 



Vn 2 3 -n? 



Fig. 2 — Asymptotic values of the surface impedances Z t and Z, for different boundary 
conditions. 
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and the coefficients X and Y are determined by the surface impedances 
given in Fig. 2. Notice there are cases where [H] is not a diagonal 
matrix, as pointed out in Section VIII. 

In the following sections we consider a waveguide with boundary 
conditions given by eq. (5) and derive simple expressions for the 
dependence of \f/(x, y) upon the matrix [H]. These expressions are 
obtained neglecting terms of order higher than 1/ka and, therefore, 
they do not require an exact knowledge of [H], but only of the 
asymptotic behavior of [H], which can be determined as seen in this 
section. In Appendix A, a procedure for determining \j/ to any desired 
accuracy is pointed out, but the procedure requires that [H] be known 
accurately. Any desired accuracy for [HI can be obtained by a proce- 
dure of successive approximations, but the resulting expressions are in 
general too complicated to be of practical interest. 

Concerning the validity of the following derivation, it is important 
to realize that even though the expressions obtained for \j/ will not 
satisfy the actual boundary conditions exactly, the errors will be small, 
of order two in 1/ka. These errors imply the conditions satisfied at the 
boundary by the expressions in question can be obtained, from the 
actual conditions, by small perturbations, of order two in 1/ka. 

III. ASYMPTOTIC PROPERTIES OF + 

We now determine the dependence of \p(x, y) upon the matrix [H]. 
To this purpose, it is convenient to assume that [if] is independent of 
ka. The expressions obtained for \(/(x, y) will depend on the coefficients 
of [H ]. By substituting for these coefficients the expressions obtained 
in Section II, one obtains the dependence of $ upon ka in general, for 
a waveguide with frequency dependent [H], 

Thus, consider a waveguide characterized by a given matrix [H], 
and let 

o 2 = k 2 -0 2 . (20) 

Let \{/(x, y) and a 2 be expanded in power series of l/ka, 

i> = M*>y)+ Z-n^M*>y)> (2i) 

«=i (ka) 

+ -*( l *k&) <22) 

where the distributions »//„, 4>i, etc., are independent of ka; they are 
determined entirely by the shape of the boundary and the coefficients 
of[H]. 

Using eqs. (21) and (22) one can derive from eq. (5) for large ka (see 
Appendix A) a set of linear relations involving \p and the normal 
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derivative ty/dv, 




4>, 




Expressing »//„ and \p r in terms of t 


fc 
fc 


6V 



+ 0[A" 2 ] on C. 



(23) 



+ 0[AT 2 ] on C, 



where 



[A]- 



— Py 



m 



— Vy 



(24) 



(25) 



v x and i^ being the direction cosines of v. 

3. 1 Derivation of i//. and C, 

We now show that for each nondegenerate eigenvalue a*, there are 
in general two modes, characterized by different values of ci . For most 
waveguides, [H] has the diagonal form (19), and in this case we shall 
see that $„, is linearly polarized. Furthermore, if the boundary has an 
axis of symmetry, then the polarization vector i of »//« is either parallel 
or orthogonal to the symmetry axis. Very simple expressions are 
obtained in this case for $* t ci,ifr.. More difficult is the treatment for 
degeneracy of order N > 1. Then, in order to determine »/>„, one must 
find the 2iV latent roots of a certain characteristic equation. 

The function ^ must satisfy the wave equation, 



Vty + <rty = 0, 



(26) 



V, being the transverse part of V. Equation (24) implies that the 
boundary values of \p vanish in the limit as k -> oo.* It is, therefore, 
convenient to represent \p in terms of the eigenfunctions f u f 2 , etc., 
that satisfy the boundary condition 



f r = on C. 
Let u r be the eigenvalue of f r , 

V?fr+U 2 rf r = 0. 

From equations (21) and (22) for k -* oo, one has 



(27) 
(28) 
(29) 



* From now on the waveguide dimensions will be kept fixed, as A' is increased. 
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and, therefore, «//„ must satisfy the wave equation with a replaced by 
a„. Furthermore, from eq. (24), 

«//» = on C. (30) 

Therefore, a x must be one of the eigenvalues u r . Let 

a. = ux (31) 

and suppose there is degeneracy of order N, so that N distinct eigen- 
functions /i , • • • , ft* correspond to the same eigenvalue. Then 

u\ = • ■ ■ = un (32) 

and \p x can be written in the form 

N N 

^ = £ CtxmfmiX, y) U + 2 OtymfmiX, y) iy, (33) 

m — 1 m — 1 

involving N eigenfunctions and 2N coefficients a xm , ciym. We now show 
that these coefficients cannot be chosen arbitrarily, but there are in 
general only 2iV possible choices corresponding to 2N distinct modes. 

3.2 Values of a xm , a ym 

The values of $ at any point (*', y') inside the boundary are related 
to the boundary values through the integral relation 24, 32 

#(*',/>- [ 4>(x,y)^dl, (34) 

Jc dp 

where G = G(x, y; x', y') is Green's function satisfying the equation 

V 2 G + a 2 G = 8(x - x')8(y - y') (35) 

and the boundary condition G = when x, y is a point of C. Let G be 
represented in terms of the eigenfunctions f r , 

G - l 2 X 2 fr(x, y)f r (x', /), (36) 

i a — u r 

where it is assumed that f r are properly normalized so that they are 
real functions satisfying 

f 2 r (x,y)dxdy=l (r=l,2, ..-), (37) 

S being the region inside the boundary C. From eq. (36), G contains a 
component 

G„ = 2 X 2 lf r (x, y)f r (x', /), (38) 

a — a oo i 
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which diverges for k — » oo, since a — > a*. For large k, 

1 1 N 
G„ = ka — -Z f r (x, y)f r (x', y'). 

Oi C] i 



(39) 



The asymptotic behavior of eq. (34) for large k is now examined. 
Approximating G by Goo, from eqs. (24), (25), and (34) one obtains the 
integral relation 



-;j>£ 



dG a 

dp 



dl, 



where 



*. 






(40) 



(41) 



Substituting eq. (39) in eq. (40), and taking into account eq. (33), one 
obtains for the coefficients a xm , a ym the characteristic equation 



Ox] 

Oty] 



Cl = 






Ox] 

«y] 



where 



and 



! Otxl 



i «yi 



(/~),,=4U"~^ 

CToo J„ W Of 



(42) 



(43) 



(44) 



and similarly for I xy , I yx , Iyy (replace Axx with Axy, A yx , Ayy). 

Equation (42) admits, in general , a total of 2N independent solu- 
tions ai], a*], -••, <X2n] for 



, OL X ] 

a] = ^ 
a y \ 



(45) 



Each solution is obtained by solving eq. (42) with C\ set equal to one 
of the 2 N latent roots Ai , A 2 , • • • , A2/V of the matrix 



[/] = 






(46) 



If the boundary is lossless, then one can verify that [/] is Hermitian 
and its latent roots are real. In this case the 2N solutions are orthog- 
onal, 



ai]«s]* = for i ^ s, 
where ( )* denotes the transpose conjugate. 



(47) 
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Thus, to determine the coefficients a x i, a y i, etc., which specify if/*,, 
the 2N latent roots of the matrix [I] must be determined. If the roots 
are all distinct, then they correspond to 2N distinct modes character- 
ized by different C\, i.e., by different propagation constant /?. If there 
is degeneracy (N > 1), the expressions obtained from eq. (33) for ^„ 
are quite complicated. Much simpler is the treatment for N = 1, since 
then only one eigenfunction f\(x, y) is involved. This is the most 
important case in practice, since it applies to the fundamental modes, 
which correspond to the lowest a m (see Appendix B). 

3.3 CaseN= 1 

For N = 1, only one eigenfunction fi(x, y) corresponds to o x and eq. 
(33) reduces to 



$*> = fi(x, v)i, 



(48) 



where i is a unit vector that determines the polarization of $„ . If Ai ^ 
X2, then eq. (42) for N = 1 specifies two polarizations, corresponding to 
two modes with different propagation constants. To determine these 
two polarizations, let oc x and <x y be the direction cosines of i. Then from 
eq. (42) with N = 1, 



where 



Ixx - Ci 



I XX = 



txy 

Iyy — C\ 



OLx 

CLy 



= o, 



dp 



dl, 



(49) 



(50) 



and similarly for I xy , I yx , ly? (replace A xx with A xy , etc.). From eq. (49) 
one obtains for C\ the characteristic equation 



(/** - C\){Iyy — Cl) = Ixylyx, 



(51) 



whose solutions Ai and A2 determine for i two eigenvectors ii and i2 
with direction cosines specified by eq. (49). Notice ii = i2 if Ai # A2. K 
the boundary is lossless and Ai # A2, then from eq. (47) 



ii-i 2 * = 0, 



(52) 



and, therefore, the two eigenvectors represent orthogonal polarizations. 
For all the waveguides of Fig. 2, [H] is given by eq. (19). Then 



Ixy = h 



r (Y~X)P X P 3 



a 

7i 



7** = -^- I (Ywl + Xvi) 



Bfx 

dp 



dp 



dl, 



dl, 



(53) 



(54) 
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and similarly for lyy (interchange x *± y). Since in this case [/] is a 
symmetric matrix, 

ii • i 2 = 0, (55) 

and if ii is real (»//„ is linearly polarized), then U is also real and 
orthogonal to ii . 

We conclude by deriving a general condition that must be satisfied 
so that \pa= is linearly polarized. Suppose i is real. Then the x axis can 
be oriented so that i = i*. This implies a y = and therefore I yx = 0. We 
conclude that linear polarization is obtained only if it is possible to 
orient the x axis so that 



I 



(Y-X)p x p > 
c 



dp 
Notice then I xy = I yx = 0, and, therefore, 



dl = 0. (56) 



ii = i*. i2 = i y . (57) 

The above requirement is always satisfied if the boundary is lossless 
(then X and Y are real) or if the values of X and Y are independent of 
position I on the boundary. It is also satisfied if the boundary has an 
axis of symmetry. In fact, let the symmetry axis be the x axis. Then X, 
Y and df\ /dp in eqs. (53) and (54) are even functions of y, 

X(x, y) - X(x, -y), Y(x, y) = Y(x, -y), (58) 

whereas v x v y is an odd function of y, which gives condition (56). 

We thus conclude that in most cases of practical interest \p x is 
linearly polarized. This is of importance in the design of a feed for 
reflector antennas, to obtain good cross-polarization discrimination 
over a wide frequency range. 

IV. PROPAGATION CONSTANT FOR N = 1 

Assume the boundary has an axis of symmetry given by the x axis 
and let N = 1. Let [H] be given by eq. (19), which applies to all 
waveguides of Fig. 2. Then, for the mode polarized in the x direction, 
the coefficient c\ coincides with I xx and it is given by eq. (54). Notice 
eq. (54) assumes that fi(x, y) is normalized as shown by eq. (37). If /i 
is not normalized, we must divide the right-hand side of eq. (54) by the 
left-hand side of eq. (37) with r = 1, then obtaining 



c, = ja — — (59) 



(Z T v l y IZ + Zpi/Z 2 ) (dfi/dp) 2 dl 
u\ [[fldxdy 
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for i = i x . For the other polarization i = i y , interchange v y *±p x in the 
above expression, which shows that the two polarizations are in general 
characterized by different propagation constants. 

Once Ci is known, the propagation constant /? can be derived using 
eqs. (20) and (22) which for a x = Mi give 

where the dots indicate terms of order higher than two in l/k. If the 
medium inside the waveguide is lossless, k is real and the attenuation 
constant -q is determined by the imaginary part of C\. Then eqs. (59) 
and (60) give 

i> = -lm(0) = ^ — , (61) 

2k 2 I I f\ dxdy 

where r and g are the real parts of Z r /Z and Z/Z z . This relation was 
used in Ref. 1 to determine the attenuation constant for a variety of 
waveguides of practical interest. 

Using the above expressions one can straightforwardly calculate the 
dispersion and attenuation characteristics of any mode for large ka. In 
the special case of a hollow waveguide of dielectric with circular 
boundary, eqs. (59) and (60) give eq. (31) of Ref. 2. 

Of greatest importance are the fundamental modes, which corre- 
spond to the lowest a* . Then, for the circular boundary of Fig. 3, 

f 1 (x,y)=AJ (o aiP ), (62) 

where p = vx 2 + y 2 , Jo is the Bessel function of order zero, and a m a is 
the first root of Jo, 

o 0B a = 2.4048. (63) 

For the rectangular boundary of Fig. 3, 

/,(,,,) = ,lcos(^)cos(^). (64) 



If fi(x, y) is normalized [see eq. (37)], then 

1 1 



A = 



ayfn Ji(<*-a) 

1 



(circle), 
(rectangle). 



(65) 



Jab 
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h— — H 






Fig. 3 — Circular and rectangular boundaries. 

V. THE DISTRIBUTION if>, 

Once i/'oo is known, the boundary values of \p can be calculated with 
error of order two in 1/k using eq. (23). If N = 1, i = i x and [H] is given 
by eq. (19), then one obtains at the boundary 



* — I 



(X v l+Yv 2 x )i x -(X-Y)p x u y i y 



dfi 

dv 



21 on C. (66) 



To determine i£ inside the boundary we separate \p into two parts, a 
component 

fi(x,y)i, 

plus a component due to the other eigenfunctions f 2 , fz, etc. The latter 
component can be determined with error of order two in 1/k by 
substituting eq. (66) in the integral of eq. (34), with G replaced by G 
— Goc for a = <7ao, as shown in Appendix A. Concerning the former 
component, it is shown in Appendix A that if the boundary has an axis 
of symmetry, then i is independent of k, and, therefore, one can set 
i = \ x for all values of k. If there is no symmetry, then in general i is a 
function of k and, to determine its dependence upon k, one must follow 
the same procedure used in this section to determine i for k—*<x>. 

VI. APPLICATIONS 

We now derive the fundamental modes of circular and rectangular 
boundaries with diagonal matrix [H ] given by eq. (19), which applies 
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to all the waveguides of Fig. 2. The surface parameters X and Y will be 
assumed to be independent of position on the boundary. Then for 
i = ix, 

^ = fi(x,y)i x , (67) 

giving «// in the limit as k — > oo. We now wish to obtain a better 
approximation for i//. 

For the rectangular boundary of Fig. 3b in eq. (66), one has v x v y = 0, 
and, therefore, i/> has the same polarization as ty m , i.e., 

neglecting terms of order two in 1/k. The boundary condition (23) can 
then be satisfied separating \p into a product of two functions, ^i(x) 
and fo(y), subject to the conditions 



Ytoh 

^! = — for x = ±a, 

k dv 

^ 2 = _ t^~ for y = ±b > 

k dv 



(68) 



(69) 



whose solutions are well known. 26 " 31 For the fundamental mode one 
obtains 



\p = A cos ax cos yy, 



(70) 



where 



IT 

aa = — 
2 



Y_ 

ka 






X 



(71) 



These results are closely related to expressions derived in Ref. 3 for a 
dielectric waveguide. Notice that for b -*■ oo the rectangular waveguide 
degenerates into two parallel plates placed at x = ±a, in which case 
yp — > >//] {x) and the modes can be derived exactly as in Appendix C. 
Similarly, for a — > oo one obtains two plates at y = ±b and ^ — » ifo(y). 
If both a and b are finite, then eq. (70) shows that ^ is simply the 
product of the two distributions ^i(x) and \fc(y), provided terms of 
order two in 1/k can be neglected. 

Next consider the circular waveguide of Fig. 3. In this case it is 
convenient to introduce polar coordinates p, <f>. Taking into account 
that X and Y are independent of <£, one obtains for the fundamental 
mode 



A 



I Jo 



(op)ix - 



X- Yja^a)' 
~~4 ka~ 



J2(ap)(cos2<|>i x + sin2<f>iy) + 



)■ 



(72) 
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where the dots indicate terms of order two in l/k, and from eqs. (22) 
and (59), 



oa = a x a 



X+ Y 1 
2 ka 



(73) 



o x a being given by eq. (63). These expressions were derived previously 
for X = in Ref. 5 and therefore details of their derivation are not 
given here. 

The above results for rectangular and circular waveguide are valid 
provided k is large enough so that the field at the boundary is small. 
From eq. (66) this requires 

1 X , Y , 

— «1, t-«1, 7-« 1 - 
ka ka ka 

If the waveguide dimensions are large enough, the results apply even 
to an ordinary waveguide with metal walls of finite conductivity. To 
determine how large the dimensions must be, consider for instance 
copper at 100 GHz. Then R s /Z = 2.189 X 10 -4 , and from Fig. 2d one 
obtains | Y \ = 514 and | X | = 5 X 10~ 5 . Therefore, the above inequalities 
require 

2a » 1000 A, 

which is too large a diameter for all practical purposes. 

The above requirement is a consequence of the large value of Y for 
copper walls. By coating the walls with a thin dielectric layer, or by 
corrugating them, much lower values of Y can be obtained. Suppose 
for instance in Fig. 2e one chooses T = 1, R 8 * 0, rti — 1, and n\ = 2. 
Then, instead of the above requirement, one obtains a »0.318X, which 
is a much more realistic condition. The attenuation constant of such 
a waveguide, or of other waveguides realized using one of the structures 
of Fig. 2, can be determined straightforwardly using eq. (61). * 

When a waveguide is used to illuminate a feed aperture, then at that 
aperture usually ka » 1. Then the aperture illumination is given 
accurately by eqs. (70) or (72) for rectangular and circular apertures. 
For other apertures, it can be determined as pointed out in Section V. 
By deriving the Fourier transform of eqs. (70) and (72) the far-field 
can be determined and thus its dependence on the aperture parameters 
X and Y can be obtained. These applications are discussed in Ref. 1. 

Helical waveguide — A case where [H] is not a diagonal matrix 

Consider one of the waveguides of Fig. 2, and let a helical wire be 
placed at the boundary, as in Fig. 4. Then one finds that [H] is not a 
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-BOUNDARY C 



K « H 

Fig. 4 — Helix with pitch angle 0. 



diagonal matrix, 

m = cos 2 0-ATsin 2 



Y, j sine eosd(l + XY) 

-j sin6 cos6 (1 + XY), X 



, (74) 



where X, Y denote the coefficients of [H] for e — > 0, e being the pitch 
angle between the wire and t. 

Consider a circular boundary. Then, we have seen that for e = the 
fundamental modes have the same propagation constant, and are 
linearly polarized for k — > <». In this section we show that for e ¥" 
they become circularly polarized, and have different propagation con- 
stants fi\ and fiz. This implies the following. Suppose at the input of 
such a waveguide (Fig. 4) the fundamental modes are combined so as 
to obtain, to a good approximation, linearly polarized excitation. Then, 
one will not obtain, in general, linear polarization at the output, unless 
the difference between the two propagation constants is small enough 
so that 

(&-/?>)/«£, 

I being the waveguide length. In this section we derive /?-2 and /?i. 
Helical waveguides are of importance for their simplicity of construc- 
tion, as compared to corrugated waveguides, and for their excellent 
performance as hybrid-mode feeds, as shown recently by R. H. Turrin 33 
whose work motivated the calculation of this section. The following 
results agree with Ref. 34. 

Let consideration be restricted to the fundamental modes of a 
lossless boundary with [H] independent of position on the boundary. 
Then, from eqs. (25) and (44) 

Lx = #1 1 V XX - {Hn + H 2i )p X y+ H-ZtVyy, (75) 

I xy = (Hu — H-22)Vxy + H\2V XX ~ H2\Vyy, (76) 
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and similarly for I yx and I yy (interchange 1 «=s 2, x +± y), where 

a/. 



Vjh = —T I VjVk 



°- , r 



dv 



dl (j,k = x,y). (77) 



Since the boundary is assumed to be lossless, 

(lm)(I xy ) = Hv2(i> xx + p yy ), (78) 

and, therefore, I xy ¥* 0. We conclude that for Hu ¥* there is no 
degeneracy possible for the fundamental modes. For a circular bound- 
ary, using eqs. (74), (25), (50), (51), and (62), we obtain 

i = -={i x ±Jh) (79) 

v/2 

and 

1 



oa = a x a 



1 - 



2(cos 2 0-XYsin 2 0) 

X+Y n „ „ 1 + XY\ 

+ 2sin0cos0 — z ! + ••• , (80) 



ka ka 

with the plus sign of eq. (79) corresponding to the minus sign of eq. 
(80). The same expressions apply also to a square aperture with a = b. 
Thus, in both cases «//„ is circularly polarized. If XY + 1 > 0, then the 
lower value of a corresponds to a mode with polarization rotating in 
the sense of the helix in Fig. 4. The opposite is true for XY + 1 < 0. 

VII. CONCLUSIONS 

To summarize, we have shown for most of the modes inside a 
cylindrical waveguide of finite surface impedances that asymptotically, 
for large values of ka, the field \p vanishes at the boundary. We have 
seen in Section III that for each eigenvalue u, there are, in general, 2N 
modes, given for k — » <» by eq. (33). For the lowest eigenvalue one has 
N = 1, and for the corresponding two modes i/^ is given by f\(x, v)i, 
where i is a unit vector. If the direction cosines of i are complex, then 
t/'oc is elliptically polarized. If [H] is a diagonal matrix, as for the 
waveguides of Fig. 2, and the boundary has an axis of symmetry, then 
i is real, and one can always orient the x axis so that i = i x . In this case 
the propagation constant (3 is given by eqs. (59) and (60), and using 
the procedure of Appendix A we can straightforwardly determine \p, 
with error of order two in 1/k, for any boundary shape. For rectangular 
and circular boundaries »// is given by eqs. (70) and (72). 

Of special importance are the fundamental modes, which correspond 
to the lowest eigenvalue o x . These modes, treated in Section 3.2, are 
needed in reflector antennas to minimize cross-polarization and edge 
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illumination over the feed aperture. They are also needed for long 
distance waveguide or fiberguide communication. Our results, show 
that there is no need to corrugate the walls of a feed in order to obtain 
conditions (11) and (48) or to obtain the low attenuations calculated in 
Refs. 20 and 21. Furthermore, they imply that the low attenuations 
predicted in Ref. 2 for a hollow waveguide of dielectric, or for a 
waveguide with metal walls coated with dielectric, 4 can be achieved 
also using other waveguides. These applications are discussed in 
Ref. 1. 



APPENDIX A 

A. 1 Derivation of eq. (23) 

Taking into account that V^E = 0, one can express E 2 in terms on 
the transverse component given by eq. (2). One obtains 

e -Jh* 
E. — V*. 

Using this relation and Maxwell's equation — junH = V X E one can 
express H 2 and H, in terms of »//. Substituting these expressions in eq. 
(5), we obtain the boundary condition 



h 



k 



1 b (V '*> 

vPJ * -V,.(i*X*)J 



1 B 



+ J^'^ 



(81) 



Taking into account eqs. (21) and (22) for k — » oo, we have 

«£->(), on C. 
This implies 

dp ov 

Taking into account these relations, from eq. (81) we obtain eq. (23) 
with error of order two in 1/k. 

A.2 Development of a 2 and ^ in Asymptotic Series of 1/ka 

A general procedure for deriving the various terms c r , «/v in eqs. (21) 
and (22) is now described, thus justifying these equations. Assume that 
the boundary has at least one line of symmetry, since this simplifies 
considerably the derivation, and it applies to most cases of practical 
interest. Also assume that [H] has the diagonal form of eq. (19), and 
that a single eigenfunction f\ (x, y) corresponds to a*,. 

Since there is no degeneracy, $„ is given by eq. (48). Let the x axis 
coincide with the symmetry line. Then, the surface parameters X and 
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Y are even functions of y, 

X(x,-y)=X(x,y) 



and 



Y(x,-y)=Y(x,y), 



(82) 



and we now show the modes can be divided in two groups, namely 
even modes satisfying 

«M*. ~y) = «M*. y), «M*. -y) = -+ y (x, y), (83) 

and odd modes satisfying 

«M«, -y) = -Mx, y), M*, -y) = «M*. y). (84) 

In fact, let 

^ = iM*, y)U + *k(x, y)\ y 

be a solution of the wave equation and of the boundary condition (81) 
with [H] given by eq. (19). Then we wish to show that 

«£' = «M*i -y)i* - »M*> -y)h 

is also a solution. Notice that \{/' is the image of »// with respect to the 
x axis, as shown in Fig. 5, where P and P' denote two corresponding 
points (x, y) and (x, —y). One can verify that 



(V,x^')p. = -(V,x^)p, 
(V^')p- = (V,^)p. 



(85) 
(86) 



If P is a boundary point, and v' and *» denote the normals to the 
boundary at P' and P, respectively, one has that v' is the image of v 
because the boundary is symmetrical. Taking all this into account one 
can verify that \f/' satisfies the boundary condition (81) with [H] given 
by eq. (19) at P'. We conclude that if an arbitrary solution \p is known, 
two independent solutions \p e and ^ can be obtained by the relations 



^e = Vfcty + 1^], +o = ViW - f ']. 



(87) 




Fig. 5 — For a symmetrical boundary, the boundary conditions are satisfied by both 
i£ and its mirror image \p'. 
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One solution, \p e , is even and the other, ^ , is odd. 

We now proceed to derive a 2 and $. Let i — i x and assume f\(x,y) is 
even, i.e., 

Mx,y)=fx(x,-y), (88) 

all other cases (i = i^, or /i odd) being entirely analogous. Then \p x and 
}p y are, respectively, even and odd and it is convenient to separate G 
into three parts, 

G = G a , + G x + G y , (89) 

where G x , G y denote, respectively, the even and odd parts of G — G».* 
Taking into account that \p x is even, from eqs. (34), (36), and (89) one 
obtains 



Jc dv 



yp x = A 1 f 1 (x,y)+ 4>x — dl, (90) 

'c 



where 



Similarly, for \f/ y , 



Ai = —z 2 r* ~T~ dl - 



(91) 



f dC 
J. = %, — & < 92 > 

Jc dv 

In eq. (90), since \p can be multiplied by an arbitrary constant, the 
coefficient Ai can be chosen arbitrarily. We choose A\ = 1, to be 
consistent with «//« = f\{x, y)i x . Then eq. (91) gives 

-ai = I + x ?p-dl t (93) 

) c dv 



~ 2 „2 _ I 

O — (Joo = I 

Jc 



a basic relation that expresses a 2 in terms of the boundary values of 
ty x . Expanding ^ina power series of 1/ka, from eq. (91) one obtains 
for the ith coefficient of a 2 , 

Ci = \\ ^«~ rf/, (94) 

°~ Jc dv 

where \p X i is the x component of ^ in eq. (21). From eq. (66), 



V.2> ^ 



^! = -(A-^-y^)a-^, (95) 

a*' 



* Thus G y is obtained from eq. (36) considering only those terms for which f r are odd, 
whereas for G It consideration is restricted to f r even with r > 1. 
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and, therefore, for i = 1, eq. (94) gives eq. (59). 

From eqs. (90) and (92), taking into account that at the boundary 

8\p = ^ — l^oc = $, 

JrtC C tlC 

Sfc—^dl, % = tyy—dl. (96) 

c dv Jc dv 

Thus, using eq. (95), we obtain 

fct - - [ (Xv 2 + Yvl) ^ &\ dl, (97) 

for the values of i// x i inside the boundary whereas for fa, using eq. 
(66), we obtain 

fa = J v x v y (x ~Y)^ (j2\ dl, (98) 

where ( )„ denotes the value for a = a*. Equations (97) and (98) allow 
«// and its derivatives to be determined with an error 0(l/k). Therefore, 
the right-hand side of eq. (81), which contains the factor 1/k, can now 
be calculated with an error 0(1/ k 2 ). The boundary values of fa, fa 
are then obtained from eq. (81). 

Once these values are known, eq. (96) can be used to determine \p X 2 
and fa inside the boundary. Equation (81) can then be used again to 
determine the boundary values of \p X 3 fa, whose values inside the 
boundary can then be calculated using eq. (96). By proceeding this 
way, we can successively calculate all fa, fa. Notice in eq. (96) that 
the kernels depend on a 2 . Therefore, to detennine fa, fa for i > 1, we 
must first calculate the coefficients ci, • • • , c,-i using eq. (94). Once 
these coefficients are known, the kernels must be developed in power 
series of 1/k, and then the first i — 1 terms in these series must be 
determined. These terms then allow eq. (96) to be used to determine 

fa, fa- 

APPENDIX B 

Let f\ be one of the eigenfunctions satisfying the boundary condition 
f x = and the wave equation 

VrYi + uif\ = 0, 

and let U\ be the lowest eigenvalue. This implies that if g(x, y) is an 
arbitrary function with continuous derivatives, then 

f f gV 2 g dxdy <«iff g 2 dxdy, (99) 
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where the inequality sign applies to any g(x, v) that is not an eigen- 
function corresponding to the eigenvalue U\. 

Condition (99) is now used to show that f\ cannot have nodal lines 
inside S. In fact, suppose f\ has a nodal line inside S, and let 

Then, since Vtg is discontinuous across the nodal line, g cannot be an 
eigenfunction, and, therefore, condition (99) should give an inequality. 
To evaluate the integral 

I=\\gV 2 gdxdy (100) 

in the immediate vicinity of the nodal line, where V 2 g diverges because 
of the discontinuity of V t g, write 

g V 2 g = V,[gV t g] - (V t g) 2 

and notice that 

is continuous because g = on the nodal line. It follows that gV 2 g 
does not diverge on the nodal line, and, therefore, its integral over a 
narrow strip containing the nodal line vanishes as the width of the 
strip goes to zero. Thus, 



1 = 



J ( fiVffi dxdy - u\ J J f\ dxdy (101) 



and, therefore, condition (99) gives an equality, which implies f\ cannot 
have nodal lines inside S. 

It is now shown that f\ is the only eigenfunction corresponding to 
U\ . In fact, if U is another eigenfunction corresponding to «i, then this 
must be true also for 

f=fi + ocf 2 , 

where a is an arbitrary constant. But this is not possible, since one can 
always choose a causing f to have a nodal line inside S, and we have 
already seen that this violates condition (99). 

APPENDIX C 

Consider the modes propagating between two parallel planes or- 
thogonal to the x axis and let 2 a be the spacing of the two planes. 
Assume at the boundary one of the conditions of Fig. 2, and let the x 
axis be oriented in the direction of propagation, so that there is no 
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variation in the y direction. Then 

$ = Hx)i, (102) 

where i is a unit vector. For the TM-modes, i = i x , whereas for the te- 
modes i = i y . For the former case E T = H 2 = and, therefore, ^ is 
independent of the transverse impedance Z T . In the latter case H T = E z 
= and \}/ is independent of Z z . In either case one finds 24,27 that the 
surface impedances Z T and Z z in eq. (9) can be determined straightfor- 
wardly. For two metal plates with dielectric coating of thickness d, for 
instance, Z T and Z z are determined entirely by /?, d, and the refractive 
index n of the dielectric, 



Y- = - fain (d Jk 2 n 2 - P 2 ) , 

Jk 2 n 2 - P 2 \ J 



(103) 



and 



y= n2k , (104) 



^^^ *n(dJZ?=W) 



which for ji -> k give eqs. (17) and (18). Analogous expressions are 
obtained for the other boundaries in Fig. 2. 
Consider the even modes with i = i x , 

$ = ipix = cosa^*, (105) 

where from eq. (81) the wavenumber a must satisfy 31 

Y 1 

(106) 



ka aa tanaa 

whose behavior for real values of aa is illustrated in Fig. 6a. If Y ^ oo, 
then Y/ka vanishes for ka^><» and, therefore, for most of the solutions 
ofeq. (106), 

aa^mir-^, m = 1, 2, • . • , (107) 

as ka -+ oo. In addition to these solutions, for Y < one mode exists 
for which aa is imaginary and 31 

oa^> -j—, as a-»oo. (108) 

Equation (107) implies that the boundary values of \p vanish for ka 
— ► oo, whereas eq. (108) implies that the mode is a surface wave whose 
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Y/ka 





Fig. 6 — Relationship between aa and Y/ka for the TM-modes of a parallel plate 
waveguide. Notice in the dielectric-lined waveguide the spacing 2a of the two plates 
increases with 2. As a consequence, the field amplitude <// at the boundary decreases 
with a, for the mode with m — 1. 



amplitude in the vicinity of the walls is given by 

k . 



^2 exp ( + ryT fl| exp 



Notice eq. (107) implies 

as ka — > 00, whereas for the surface wave 



x — a 



fi-+k 



> + ^*. 



(109) 



(110) 



(111) 



To understand the significance of these results, consider two metal 
plates with dielectric coating as in Fig. la. Let the separation a of the 
two plates be gradually increased in the direction of propagation, as 
shown in Fig. 6, and let a — > 00 as 2 — > 00. Let the dielectric thickness 
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d be so small that initially, for 2 = 0, 

— = - 00. (112) 

ka 

Then the modes for z = are essentially those of an ordinary waveguide 
without dielectric coating, 

oa=*mir (m = 0, 1, 2, • ••)• (113) 

For z > 0, however, the magnitude of Y/ka decreases with z, and 
Y/ka — > for z -* oo. This implies, for all the modes with m ^ 0, 
conditions (107) and (110) and, therefore, the boundary values of i// 
vanish for z — » oo. For m = 0, on the other hand, aa is imaginary and 
initially aa = 0. This mode will degenerate for z — » oo into a surface 
wave with propagation constant determined by Y as shown by eq. 
(111). This is the only mode for which the field does not become 
infinitesimal at the walls for z -» oo. For all the other modes the 
boundary values of i// for large ka are given by 

(-ir +1 Y^, (114) 

ka 

which vanishes for ka — > oo. 

The above considerations apply also to the TE-modes. In fact, it can 
be verified that if in eq. (2) E t is replaced with ZH, so that \p represents 
the transverse component of ZH t , then, for the even modes, \p is still 
given by eqs. (105) and (106), provided Y is replaced with X. The 
behavior of the odd modes is entirely analogous; simply replace 
cos ax with sin ax in eq. (105), and tan with - cotan in eq. (106). 
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